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Abstract 

According to modern developments in turbulence theory, the "dissipation" scales (u.v. 
cut-offs) r\ form a random field related to velocity increments (Lit. In this work we, 
using Merlin's transform combined with the Gaussain large -scale boundary condition, 
calculate probability densities (PDFs) of velocity increments P(5 r u, r) and the PDF of 
the dissipation scales Q(rj,Re), where Re is the large-scale Reynolds number. The re- 
sulting expressions strongly deviate from the Log-normal PDF PL(5 r u,r) often quoted 
in the literature. It is shown that the probability density of the small-scale velocity 
fluctuations includes information about the large (integral) scale dynamics which is 
responsible for deviation of P(5 r u,r) from Pi,(S r u,r). A framework for evaluation of 
the PDFs of various turbulence characteristics involving spatial derivatives is devel- 
oped. The exact relation, free of spurious Logarithms recently discussed in Frisch et 
al (J. Fluid Mech. 542, 97 (2005)), for the multifractal probability density of velocity 
increments, not based on the steepest descent evaluation of the integrals is obtained 
and the calculated function D(h) is close to experimental data. A novel derivation 
(Polyakov, 2005), of a well-known result of the multi- fractal theory [Frisch, "Turbu- 
lence. Legacy of A.N.Kolmogorov v , Cambridge University Press, 1995)) , based on the 
concepts described in this paper, is also presented. 
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1 Introduction. 



A reasonably well experimentally established anomalous (multi) scaling of the structure 
functions S n>m = (u(x + r) — u(x)) n (v(x + r) — v(x)) n = (5 r u) m (5 r v) m , where u and v are 
components of the velocity field parallel and perpendicular to the displacement vector r = ri, 
respectively, is one of the properties of strong turbulence which makes it "the last unsolved 
problem of continuum mechanics". [l]-[4]. The anomalous dimension, a property of strongly 
interacting systems first introduced to the quantum field theory by Gribov and Migdal [5] and 
Polyakov [6], is a notoriously difficult from theoretical viewpoint concept and only recently, 
after many years of trying, the theory of anomalous scaling of the structure functions of a 
passive scalar advected by a white- in-time random velocity field, has been developed [7]- [8]. 
Almost simultaneously, the theory of bi-scaling in turbulence generated by the forced Burgers 
equation [9]- [10] was formulated. The attempts to explain anomalous dimensions in three 
dimensional turbulence were made in Refs.[ll]-[12]. It was shown that the Navier-Stokes 
equations combined with a simple model for the pressure -velocity correlations , lead to 
homogeneous differential equations for the structure functions and, as a result, to anomalous 
scaling exponents which cannot be obtained on dimensional grounds. 

During last forty-fifty years, a substantial effort was devoted to derivation or at least modeling 
of various probability densities in turbulence. The first attempts resulting in the Log-normal 
PDF of the dissipation rate fluctuations, consistent with the multi-scaling, were made by 
Kolmogorov [13] and Yaglom [14], using a simple cascade model. This result was later 
criticized by Orszag [15] as, in general, not realizable. Similar Log-normal PDF was obtained 
for the not too large magnitudes of velocity increment 5 r u = u(x + r) — u(x) in Ref. [11] , 
which was later experimentally tested by Kurien and Sreenivasan [3]. Some other attempts 
based on analysis of experimental data led to various fits ranging from Log-normal and Log- 
Poison expressions to exponential and "stretched exponential " PDFs P(S r u). In this paper 
we, addressing this problem, will restrict ourselves by considerations based on equations of 
motion. 

This paper is organized in a following way. In the Section 2 we introduce the Mellin trans- 



2 



form for the probability density of velocity increments and define the large- scale Gaussian 
boundary condition. We show that the " normal (linear ) scaling" corresponds to the Gaussa- 
ian PDF of the small-scale fluctuations. The PDF accounting for small deviations from the 
linear scaling is calculated and compared with traditional Log-normal expressions. In Section 
3, we based on theoretical analysis of the Navier-Stokes equations and numerical simulations 
by Gotoh and Nakano [16], justify the boundary conditions used for the derivation of Sec 2. 
The random field of the dissipation scales rj (linear dimensions of the dissipative structures) , 
derived from the expression for the dissipation anomaly, is introduced in Sec 4 and the PDF 
Q(i], Re) is computed in Section 5. Discussions and conclusions are presented in Section 
6. There, the exact expression for the probability density of velocity increments following 
the multifractal formalism is presented and one of the results of the multi-fractal theory is 
derived (Polyakov 2005) without the multifractal input. Some of the concepts and relations 
presented below have been reported in the previous publications [12], [17]- [18]. Wherever 
needed, we repeat them here for the sake of continuity and clarity. 

2 The Mellin transform. 

If the moments of velocity increments 5 r u are given by the scaling relations S n fl = (5 r u) n = 
A(n)r^ n , then probability density function can be found from the Mellin transform: 



where we set the integral scale L and the dissipation rate S equal to unity. Multiplying (1) 
by (5 r u) k and evaluating a simple integral, gives S^o = A(k)r^ k . 

To make use of the relation (1) the dynamic information about both the amplitudes A(n) 
and the exponents £ n is needed. In order to obtain an expression for A(n) from a large 
scale boundary condition on the PDF, we first have to define the integral scale L. Based on 
experimental data and theoretical consideration (see below), we choose L as a scale at which 
the energy flux toward small scales changes sign or tends to zero. This means that at the 
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small scales r < L the structure function S3 (r) < and at the larger scales r > L, S3 > 0. 
Typically, at this scale which depends upon the geometric details of the flow, the odd-order 
moments S2 n +i(L) = and the even-order ones saturate: d r S2 n ,o(L) = 0. This property 
of turbulence has recently been examined both numerically and experimentally [19] -[20] (a 
theoretical argument will be given below). The scale L appears naturally in a simple case 
of the Navier-Stokes equations on an infinite domain driven by the white-in-time forcing 
function with the variance: 

7^= 5|^<* -*/)/*" 

where V is the forcing power. The exact calculation of the relation for the third-order 
structure function S 3 (r) gives an oscillating expression 

— 36r cos r + 12 sinr — 12(— 2 + r 2 ) sinr 
S3 = -V 

r 4 

In the limit r — > 0, we have the Kolomogorov relation S3 = —-^Pr. At the large scale 
L ~ 5.88/ kf , S 3 (L) = [21]. In all flows studied, the probability density P(5lu,L) was 
extremely close to the Gaussian. We would like to stress that the integral scale defined 
this way is not the largest scale (size) a system but rather corresponds to the top of the 
inertial range where the constant energy flux toward small scales sets up. The Gaussian 
boundary condition at r = 1 leads to a plausible and well- tested (both experimentally and 
numerically), expression A(n) = (2n — 1)!! which will be used in all calculations of this paper. 
In what follows a dynamic and numerical justification for this result will be presented. 
First we consider the case of "normal scaling " £ = an. Writing (2n— 1)!! = e~ x2 x 2n dx 
and rotating the integration axis by 90°, we have setting 5 r u = u: 



P(u,r) = i°° e' x2 dx [°° e m( lnr ^+ ln x )dn = [°° e~ x ' ''5( lnO^- + lnx)dx (2) 

^TTU J -00 J-00 V 7 ™ U 

This integral is evaluated readily with the result: 

1 u 2 

P(u) = —^= — e^zT^) 
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Now we consider anomalous scaling by introducing small deviations from the linear relation: 

= an — bn 2 (3) 



which for not too large moment numbers n can perceived as first two terms of the Taylor 
expansion of £ n in the vicinity of n = 0. The formula (3) is a generic perturbative expression 
for the exponents £ n independent upon the nature of the problem and the relations similar 
to (3) have resulted from the recent perturbative theories of a passive scalar in a random 
velocity field. [7] -[8]. Using the Kolmogorov constraint £ 3 = 1, gives: b = (3a — l)/9. It is 
clear that the expression (3) cannot be correct for all values of the moment order n. Indeed, 
in accord with Holder's inequality, £ n , is a concave and non-decreasing function of n or in 
other words as n — > oo, £(n)/n — > (£ n /n ^> 1/n). Still for n < 10 — 15, the experimental 
data on strong turbulence are consistent with a ps 0.383 and b 0.0166 and the expressions 
derived below can be accurate only for not-too-large values of velocity increment 5 r u. Thus, 
the probability density is given by the integral: 

P(u,r) = -L- r e~ x2 I" e in in^-bnHn r dn (4) 

y/TTU J-oo J-oo 

which is reduced to: 

The integration over n leading from (4) to (5) was based on the following estimate. Since 
e %n\n{u/r a ) j g an oscillating function, the main contribution to the integral over n comes from 
the interval n ~ ^/ \Jb\ \nr\ and in order for the relation (3) to be valid in the integration 
interval, the following condition must be satisfied: 1/ \Jb\ \nr\ < 10 — 15 which is satisfied 
when the displacement r is small enough. 

The PDF P(^,r) numerically evaluated from equation (5), is shown on Fig. I. for a few 
values of displacement r. We can see that the tails of the PDF strongly depend upon r which 
is a sign of intermittency and anomalous scaling. 

In the range u/r a p^ 1, the integral is dominated by the interval x ~ 1 leading to the 
log- normal result: 

5 



Both PDFs (5) and (6) are normalized to two. The expression (6) , which has also been 
derived directly from the Navier-Stokes equations combined with a simple model for the 
pressure -velocity correlation function in Ref. [11], has been experimentally verified by 
Kurien and Sreenivasan [3]. The comparison between (5) and (6) showing a surprisingly large 
difference between the approximately and numerically calculated integral (5) is presented on 
Fig. 2. In addition, it follows from Figs. 1 and 2 that while P(5 r u,r) has a maximum 
at 5 r u = 0, in the limit 5 r u — > 0, the Log-normal PDF Pl, given by (6) rapidly decreases 
to zero, meaning that it may be a reasonable approximation not too close to the origin at 
5 r u = 0. 

In the most interesting case = j^-, the integral (1) with A(n) = (2n — 1)!! can be 
evaluated both exactly and by the steepest descent method leading at the large magnitudes 
of the argument u/r a to the algebraically decreasing probability density consistent with the 
saturation of the exponents £ n — > a/j3. 



3 The relations between moments. Simulations by Go- 
toh and Nakano (Ref. [16].) 



In this Section we would like to develop some dynamic arguments justifying the choice of 
the structure functions amplitudes A2 n ,o = (2n — 1)!! used in the above calculations. The 
relations for the moments of velocity difference S n>m were derived in Refs. [11], [12] and later 
in [22] using an alternative approach. In particular, the equation for the even-order moments 
is: 




+ 



d-l 



(d- l)(2n- 1) 



S 2n -2,2 = (2n - 1) < S r au n 2 > 



(7) 



dr 



r 



r 



and taking into account that the dissipation contribution can be neglected [11], [12], [3], [22]: 
we have 




+ 



d- 1 



(d- l)(2ra- 1) 



5 , 2n-2,2 = -(2n - 1) < S r (d x p)u 



2n-2 



> 



(8) 



<9r 



r 



2n,0 ~ 



r 



This relation is not closed. In the vicinity of the integral scale where the odd-order moments 
S 2n +i(L) and derivatives of the even order moments d r S2 n ,o(L) disappear, the pressure gra- 
dient terms can be neglected. Here we would like to present a possible physical explanation 
of this effect. If geometry of the large- scale unstable structures leading to the energy cas- 
cade toward small scales is defined by the condition S 3 (L(x,y, z)) = 0, then associating the 
instability with a flow separation phenomenon, gives the standard relation (see any textbook 
on hydrodynamics): Vp = uV 2 u, meaning that at these points a = and 5l& = 0. This 
mechanism, if locally correct allows one to neglect acceleration contributions to the equation 
(7) (or pressure gradients in (8)) with the remaining terms giving S^o — (2?i — l)<S2n-2,2j 
consistent with the Gaussian distribution. Thus, extrapolating this result into the "inertial 
range", we assume that (2n — l)S 2 n-2,2 = S 2n ,oSo,2/S2,o — (1 + y^n.o which is consistent 
with the Gaussian distribution at r = 1. This relation can be proven as follows: since uv = 0, 
at a gaussian point r = 1 we have S 2n -2,2 = S , 2 S 2n -. 2>0 = ^pi^^n.o = ^-f ^n.o- One 
remark is in order: in the inertial range £2 7^ while in the limit r — > L, the r-derivatives 
disappear or in other words loosley defined in this limit "exponent" £2 — ► 0. Thus the above 
relation has to be treated as a parametrization valid in both limits r — > L and r/L <^1. In 
their remarkable paper Gotoh and Nakano [16] made a detailed examination of the various 
terms in the equations (7) and (8). Some of their results are presented on Fig. 3. where the 
left side of equations (7)-(8) denoted as "lhs" is compared with the right side ("rhs" minus 
pressure-gradient terms) for n = 4; 6; 8. On the Fig. 3 A stands for the Taylor micro-scale. 
They demonstrated that : 1. indeed , the dissipation contributions to the equations for 
the even order moments are negligibly small; 2. the relation (2n — 1)5 , 2 „_2,2 ~ 1.35<S2 n ,o, 
independently upon the displacement magnitude r (not shown on Fig. 3); 3. The pressure 
contributions (2n — l)5 r (d x p)(5 r u) 2n ~ 2 / lhs ~ const in the entire range of the displacement r- 
variation. Substituting (2n— l)S , 2n _2,2 ~ l-35<S2n,o an d the algebraic relation S^o — ^2n,o r ^ 2n 
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into ( 8) gives (d — 3): 



(dr + g)g2n,0 - frn ~ 6 (2fl - 1 ) < ^ > 

and for £ 2 ~ 0.7, £ 4 fx 1.27, £ 6 fx 1.77 and £ 8 fx 2.15 the magnitudes of the ratio: R 4 = 
0.18; R 6 = 0.28 and i? 8 = 0.35 in a close agreement with numerical data of Ref. [16] (For 
comparison, see Fig. 3). Moreover, one can see form Fig. 3 that all pressure-gradient velocity 
correlations become negligibly small in the vicinity of a single length scale (L fx lOOOr/^) 
consistent with the above introduced definition of integral scale L. Here r] K is the Kolmogorov 
scale. The equations (8) combined with the above considerations yield another interesting 
relation for the pressure gradient-velocity correlations: 

(6n - 6) Vo = "(2n - l)r5 r (d x p)(5 r uy n - 2 > 
The most important outcome of the formula (9) and numerical data of Ref. [16] is that the 
ratios i?2n = const are independent upon the magnitude of displacement r. This means that 
the Lagrangian acceleration contribution to (8) can be expressed as a linear combination: 

— < 5 r d x p\5u, 5v >= a p d r (S r u) 2 + P p (S r u) 2 /r + 7 p (5 r f ) 2 /r + K p d r (5 r v ) 2 (10) 
plus sub- leading terms which are irrelevant in the small-scale limit r — > 0. 



4 Dissipation anomaly. Calculation of spatial deriva- 
tives 

Similar result is obtained if we consider the dissipation anomaly introduced by Polyakov [9] 
who was interested in turbulence generated by the random-force-driven Burgers equation. 
Later Polyakov's results have been generalized to the Navier-Stokes equations by Duchon and 
Robert [24], Eyink [24] and Yakhot and Sreenivasan [18]. Interested in various characteristics 
of the velocity field, one often needs to calculate spatial derivatives defined as usual: 
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u(a; + y) - u(a;) 



2/^0 y 

The problem is that if the velocity field is differentiable, then S3 = 0(y 3 ) and d y S 3 — > 0. 
This is in contradiction with the Kolmogorov relation stating d y Ss oc £ = 0(1). This means 
that in the "inertial range" the velocity field is singular. To resolve this problem, Polyakov 
[9] introduced a limit-ordering procedure: 1. v — > 0; 2. y — > 0. Assuming that there exist a 
scale (" dissipation" ) separating analytic and singular ranges, we can approximately redefine 
derivative as [18]: 

d x u(x) — lim (u(x + y) — u(x — y))/2y ~ + 77) — u(x — i]))/2r] (11) 

This definition leads to an interesting consequence. In the analytic interval u(x + r/) pa 
u(x — 77) + 2i]d x u the matching condition at the scale 77 reads: ~ 2i]d x u thus defining 77 
as a random field. Therefore, to calculate the moments of derivatives, one needs either the 
probability density P(S v u, 77) or the relation coupling 77 and 5 v u. The dissipation anomaly 
enables one to calculate this relation. Here we present the main steps of derivation of Refs. 
[9], [24], [18]. 

We are interested in the Navier-Stokes dynamics of incompressible fluids, for which the 
energy balance equation (with the density p = 1) is written as 

1 du 2 1 _ 9 _ <9 2 u 
2^ + 2 U ' VM =- Vp - U + UU -dxJ> 

The differential equation for the scalar product u(x + |) ■ u(x — |) = u(+) ■ u(— ) can 
be written as 

^±^ + u(+) .^ u(+) . u( _ )+uH .^ u „. u(+) . 

-§£f"'<-> - 1£r"' (+) + ^ |u( - ) ■ aS- u(+) + u(+) ■ sb uH| - (12) 

It is clear that in the limit y — > 0, for which x ± — > x, this equation gives the energy balance. 
In the limit y — > 0, the equation (12) has two kinds of terms: the regular ones which 
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disappear by the virtue of the energy balance and a few singular terms balancing each other. 
The calculation presented in detail in Ref. [18] leads to the equation (See also Duchon and 
Robert [23], Eyink [24]): 

d Id d 

lim y ^ [-—5u l \5 y u\ 2 + 2(-^— u *(+) u j(-) 2 + fa~ Ui ^ u ^ + ^ = ~ 25 y u ' ^ ( 13 ) 

where a = — Vp + z/V 2 u is the lagrangian acceleration. The above equation, not involving 
time derivatives, is exact at each point in the flow at each instant of time. Choosing the 
displacement vector along one of the coordinate axes and averaging, one obtains [24]: 



8 d 2 

— Su\5u\ 2 = 8u5ui— 
ay ay 



5u\5u\ 2 = 8u5ui—5ui = 2iy(S y Ui)d 2 (5yUi) = --£, 



where 5 y u = 5 y u ■ y/y. The pressure terms in and the second contribution disappeared by 
the averaging procedure. 

Dissipation anomaly as a closure. The relations (7)-(8) are valid for all magnitudes of 
the displacement vector r including r — > rj. Substituting the expression for the dissipation 
anomaly (13) into the right-side of (7) gives: 



^ + ^s 2 „ - {d - l)i2 "- l) s 2n ^ = -(*, - 1) < mv)u«'--> >■ 

or r r 



3(2n - 1) < (gp» + _ d^vf > 
4 or or 

The velocity-kinetic energy product, which appears in equation (13) has not been included 
in this relation as small. Using the gaussian boundary conditions implying (2n — l)5 , 2 n _2,2 — 
(1 + ^2/2)5 , 2n,o this expression can be somewhat simplified. The full closure is achieved if we 
express the second contribution to the right side of (14) in terms of S r u. At this point we 
cannot do it rigorously but it is plausable to assume that it is 0((5 r u) 3 ) leading to a generic 
relation (10). A simple closure based on the Bernoulli-like relation 5 r (d x p) = 0(d r (5u) 2 ) 
leading to the one-parametrical expression for the exponents £ n in a close agreement with 
available experimental data has been proposed in Refs. [12]., [16]. 
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5 Dissipation structures. PDF of dissipation scales. 



On the dissipation scale rj all contributions from pressure, advection and dissipation terms 
are of the same order and the relation for the dissipation anomaly allows the estimate [18]: 

i]5 v u fx v (15) 

which means that the dissipative structures correspond to the local magnitude of the Reynolds 
number 

Re v = r]5 v u/v = 0(1) (16) 

This defines rj as a random field. We would like to stress that the relation (15) has been ob- 
tained by balancing various contributions to the locally exact expression for the dissipation 
anomaly and by establishing the relation between the dissipation scale and velocity incre- 
ment, this formula enables one to evaluate spatial derivatives and compute various correlation 
functions. If the displacement y is in the analytic range, then u(x + y) — u(x) fx d x u(x)y and 
extrapolating y — > rj where rj is the scale where the analytic and singular ranges overlap, we 
obtain [18]: 

d x u fx S v u/i] fx (5 v u) 2 jv (17) 

acceleration 



5 v u _ (5 v u) 2 _ (5 v uf ^ 



T Tj V 



and the dissipation rate 



£ = u5 v uV 2 5 v u fx (5 v u) 4 /u (19) 

It follows from (17) (See also Ref. [18]) that for each moment S n ^ m one can define a "dissi- 
pation scale r] n+m separating analytic and singular ranges and 



7] n fx LRe^-^+i- 1 (20) 
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Using the relations (17)- (20), we can develop the multi-scaling algebra. For example, 

tfH » (-^f-S 6n (rj 6n ) oc ( -) 2n rjt ~ ( —y^- ) 2n Re a2n , (21) 

with a 2n = 2n + 7 — I 52 7. With £6 = 2 and £7 = 7/3, we recover Yaglom's result [241 

9 
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a 2 ~ ^^jr- If evaluating the exponents, one uses in (21) the anomalous exponents ^% ~ 1-77 
and £7 ~ 1.99 from Ref. [12], then a 2 ~ 0.55 meaning that the intermittency correction is 
ps 0.05. Recent experiments by Reynolds et al. [25] were in a close agreement with (21). A 
similar formula has recently been obtained from the multi-fractal formalism by Biferale et al. 
[26]. Formula (21) shows that the second moment of Lagrangian acceleration is expressed in 
terms of the sixth-order structure function evaluated on its dissipation scale To extract 
information about the fourth moment a 4 , we have to have accurate data on 612(7712) which, 
in high-Reynolds-number flows, is very difficult to obtain from both physical and, especially, 
numerical experiments.. 

The moments of velocity derivatives evaluated easily. In accord with (17) and (20): 



(d x u) 2n ~ (^) 2n PS (^y-) 2n « Re 2n S± n (r) 4n ) ps Re d2n , (22) 



J4 



where d 2n = 2n + , , . . 
Moments of the dissipation rate. The exponent fi 2 . Using (19) and (20), the moments 
of the dissipation rate are calculated readily: £ n ps Re^ n where \i n = n + £4 n /(£ 4ri — £ 4n+ i — 1). 
Taking in accord wit Refs. [[11], [12] £ n ps 0.383ra/(l + 0.05n), gives 

fx 2 ~ 0.16 

and, since according to (20) i] 8 ps i?e~ 84 , we derive S 2 ps % °' 19 . This result is based on the 
relation (19) telling us that the mean of the square dissipation rate must be evaluated in terms 
of finite difference (velocity increment) over the region of space having the linear dimension 
r] 8 which is substantially smaller that the Kolmogorov scale 0(i?e -3 / 4 ). Extrapolating this 
relation into an inertial range gives 



£{x)£{x + r) oc r M 
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with /i 0.19. The same expression, calculated on the Kolmogorov scale gives /i ~ 0.21. 
It is interesting that the exponents ji n for n < 1 are negative. This prediction is yet to 
be tested experimentally. Evaluation of the high-order moments of the dissipation rate 
involving the correlation functions S in {rj An ) in their respective analytic intervals requires 
very high resolution of the velocity field and is highly problematic. 

Dissipation scales. According to its definition, the dissipation scale is a linear dimension 
of a structure defined by the local value of the Reynolds number Re v = r/S^u/u = 0(1). 
This introduces a random field. The probability density Q(i], Re, Re v ) is found by fixing 
the displacement r = i] and counting the events with r\b r u/v = Re v keeping the global 
and local Reynolds numbers Re = u rms L/u and Re v fixed as parameters. Defining r\ as 
a scale at which the analytic and singular parts of velocity field overlap, naturally leads 
to the matching condition: P(5 v u/r] a ,ri) = P D (5 v u,r]) where is the PDF of 5 r u for the 
displacement values from the "dissipation" range r < rj where the expression (3) does not 
work. On the matching scale however, the relation (3) is correct. Thus, the probability 
density Q(r),Re) can be found from P{5 7) u\Re r) m 1) where Re v = "qb^u/v. From formula 

(1): 

P(5 v u) = P(u v ) = — / A(n)u^ n) u-^ n) ~ n dn (23) 

and for the scaling exponents given by (3), the result is derived readily. Introducing the 
large-scale Reynolds number v = 1 /Re and taking into account that Reu v ~ 1 / 77 gives for 
the probability density Q{rj) : 



Q(r),Re) = - f° e^ dx f°° ^™M^ +1 ^e)-&n 2 inr, (24) 

Tj J— 00 J— 00 



l r°° 2 c°° 
Q{V,Re)= ' ""' 

giving: 



Q(rj,Re) = — . = , / e x dxe 4 « n " (25) 
•q^Abmr] J-00 

The probability density Q(rj, Re) is plotted on Fig.4a for Re — \/v — 0.3 and the curves for 
Re = 0.1; 0.3 are compared on Fig. 4b. 
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PDF of velocity derivative is computed from (23) combined with (17). The result is: 



The integral is evaluated using the procedure of Section 2. It is interesting that for a "normal" 
case £ n = n/3 we have: 



This relation, which has been obtained by Benzi et al [27] from the multi-fractal formalism, 
cannot be a consequence of the Kolmogorov (K41) theory. Indeed, unlike the expression (17) 
based on the idea of the fluctuating u.v. cut-off, the K41 dissipation scale t)k = const is a 
number not related to 5 v u . Thus, since the fluctuating dissipation scale is not compatible 
with K41 and "normal scaling", this PDF is not realizable. If for consistency with K41, one 
uses the constant dissipation scale, then the PDF is simply the one of 5 m u which, in accord 
with calculation of Sec 1, is a Gaussian. 

The PDF of the dissipation rate can be easily calculated from the integral (23) with u v xs 
(£z/)i. The resulting expression gives P e (£,Re) with the broader tails than those of Kol- 
mogorov' s Log-normal probability density. 

6 Discussion and conclusions. 

The theory presented in this paper is based on two principle assumptions : 1. The existence 
of a length-scale L such that for r = L the odd-order moments S 2n +i{L) ~ 0. At this 
scale the energy flux toward small scales sets up and the PDF P(S r u,r = L) is close to the 
Gaussian. This statement , consistent with the Navier-Stokes equations, has been tested 
in both physical and numerical experiments [19]- [20]. The Gaussian PDF P(5 L u,L) is not 
unlike a large-scale boundary condition needed to solve the differential equations (7)- (8) for 
the moments S njjn . 2. The expression (3) for the first two terms of the Taylor expansion 
of the function £„, consistent with the Holder inequality, is a good approximation to the 
exponents of the first 10-15 moments of velocity increments . The existence of a small 




P(u') oc e 
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parameter in turbulence theory is highly problematic and numerical smallness of deviations 
from the linear expression for not too large moment numbers n in (3) was helpful for the 
theory developed above. The relations for the dissipation structures and " dissipation" scales 
i] coming from the order-of-magnitude balancing of the terms in the exact equations for the 
dissipation anomaly are well justified. 

Using these assumptions we have shown that the calculated PDF of velocity differences 
strongly deviates from the Log- normal distribution, first obtained by Kolmogorov in 1962, 
which even today is widely used in the literature. This difference stems from contributions 
of the amplitudes A(n) to the integral (1). The amplitudes A(n) are fixed by the large-scale 
boundary conditions, leading to a natural conclusion: the small-scale dynamics are strongly 
coupled to the large-scale phenomena. This may be a reason for a serious reexamination of 
the very concept of the turbulence energy cascade which, within the framework of the present 
development, seem neither possible nor needed. An accurate experimental and numerical 
comparison of the measured and Log-normal PDFs of velocity increments may be extremely 
important. 

It has been shown [18] ( in a different way this is also an element of the multi-fractal theory) 
that the scales rj form a random field not necessarily related to the energy dissipation scales 
but rather to the linear dimensions of various dissipation structures defined by the local 
value of the Reynolds number Re v = 0(1). Some of these structures are responsible for 
the energy and the second-order moment S 2t o dissipation, while others, more powerful, for 
the dynamics of the higher -order moments. Thus, the scale r] must be perceived as a 
dynamic cut-off separating analytic and singular components of the velocity field. The 
probability density Q(i],Re), calculated in Section 5 is an interesting and easily measurable 
quantity. A note of caution is in order: to make reliable calculations or measurements of the 
moments involving spatial derivatives of velocity field, one has to have a field resolved well 
enough to exhibit at least a fraction of the analytic range of the corresponding moments of 
velocity increment. For example, according to (18), the second order moment of Lagrangian 
acceleration is proportional to the sixth order structure function calculated on the scale r) e , 
while the fourth order one is expressed through Su (^12)- In the high Reynolds number flows, 
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the measurements of the twelveth order structure function including analytic range where 
oc r 12 do not exist. The situation with the moments of dissipation rate is even worse: 
fourth- order moment £ A is related to Si 6 (?7i 6 ). An interesting possibility is being explored 
by Schumacher [28 ] running very large (1024 3 ) numerical simulations at reasonably low 
Reynolds numbers R\ m 10 — 60. Analyzing the probability density of the dissipation scales 
Schumacher and Sreenivasan [29] obtained Q(rj, Re) very similar to one shown on Fig. 3. 
The results presented here were obtained from analysis, both theoretical and numerical, of 
the dynamic equations. No multi-fractal assumptions have been made. Still, it is interesting 
to compare the two approaches. In its present form, the multi-fractal (MF) theory consists 
of two parts [1]. The first one based on an idea of fractal dimension, attempts to explain the 
origin of anomalous scaling by assuming 

S p (r) = = / ^)(^r +3 - D(/l) (26) 

The normalised structure functions S p = (^-) p - The 0((r/L) ph ) term comes from the 
multi-fractal assumption 

5 r u = (r/L) h (27) 

defined on a set of fractal dimension D(h) where h is a value of the scaling exponent from 
the interval h min > h > h max and (r / L) 3 ~ D ^ is the probability of being in the interval r in 
a volume of dimension 3 — D(h). Neglecting the Logs in the steepest descent evaluation of 
the integral (26), gives a relation between the fractal dimension and the exponents (For 
the review see Frisch [1]). 

The multifractal PDF and Mellin transform. Let us establish possible relations be- 
tween the two theories. Using the gaussian expression for the amplitudes, we have: 

S 2p (r) = (2p - = (2p - 1)!! / d^h)^)^' ^ (28) 

Substituting this into (1) and repeating the calculations of the section 2 gives (L — 1): 

P(u,r) = I [e-*dx r dn [ iM dhe ^it-^%^W)^ 
v ' ; 2 J J-ioc J dh 
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Integration over n gives the delta -function 5(h\nr — (u/x oc r h ) with the finial 

result: 

1 /" j__-x a ___/o r^ ln v^n^^*) 



In " 

where /i* = . Restricitng ourselves by the relation (3) for the exponents and comparing 
this formula with the PDF (5) gives (neglecting the //-factor ): 

>77§7^ (In: 



[3 - D(-^f )] lnr = - ^f* (29) 
| lnr I 46| lnr| 

No steepest descent approximation has been used in deriving this relation. If in accord with 

the MF theory, we set the amplitudes A{n) = 1 and \/2x = 1 in (29) and taking into account 

that h = In uj lnr, the expression (29) gives: 

6 UW ~ 46(lnr)2 " 46 
We remind the reader that in accord with Ref. [11], [12], a 0.383 and b ss 0.0166. This 
derivation which did not involve the steepest descent evaluation of the integral does not 
have a "Log-problem", discussed in Ref. [30]. The experimental measurements of Cramer 
function f(a) = D{h) + 2 for h = a/3 by Meneveau and Sreenivasan [31] (See also Ref.[l]) 
are in an extremely close agreement with this expression. The quantitative differences are: 
the maximum of the calculated curve (with a = 3h) is at a — 1.15 (instead of a — 1 of Ref. 
[31]) and f(a) — at cti — 0.369 and a 2 = 1.92 compared with a,\ 0.5 and a 2 ~ 1.8 of 
Ref. [31]. The small deviations come from the difference between h = a/3 used for analysis 
of experimental data and the theoretically obtained h = a = 0.383. This difference decreases 
if the coordinates are rescaled by factor 0.383/0.333. 

The second part of the MF theory, dealing with the small-scale properties of turbulence, is 
based on the relation (Paladin and Vulpiani [32] ) 

y « Re~^K (30) 

obtained by combining the MF assumption (27) and the outcome of the balance of the 
advective and viscous contributions to the Navier-Stokes equations. All small-scale results 
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derived using the MF formalism are numerically indistinguishable from the ones obtained 
both above and in the Ref. [18]. To illustrate this point, we present an alternative derivation 
of the moments of velocity derivative due to Polyakov [33] . 

Polyakov's derivation. The probability density of velocity difference in the inertial range 
is: 

P(5 r u, r) =< 5(5 r u - [u(x + r) - u(x)\) > (31) 
and in the dissipation (analytic) range 

P D (5 r u, r) =< 5(5 r u - u'r) > (32) 

Assuming the two PDFs match at the scale r — rj — v/5 v u = u/u, introduced by (15), we 
have: 

P D (u,r]) = P(u,rj) =< S(u - u'{x)v/u) >= J dnA{n)v^ n u'^ n ' 1 (33) 
Multiplying (33) by u 2k and integrating over u gives: 

I r- i 

-(yu') k = dnA(n)v in — 

The integral is evaluated at a pole where + n = 2k giving 

{u'{x)) k oc A{n{k))v p W (35) 

with 

P(k) = Cn(fe) —k = k — n(k) (36) 

which is identical to the formula (8.76) of Frisch's book [1] obtained using multi- fractal 
theory. Comparing the relations (35) and (22) we find that , on the accepted magnitudes of 
exponents £ n , numerically they are basically identical for not too large moment numbers n. 
It is also easy to see that if in the limit n — > oo, £ n oc n a with < a < 1, the two relations 
have the same asymptotics. 
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(34) 



The only approximation involved in derivations of both relations (22) and (35), presented in 
this paper, is the choice of the cut-off i] = u/u instead of rj = 0{v/u). In reality, there exist 
a random field of the dissipation scales described by the probability density Q(r),Re) given 
by (25). Thus, a more accurate calculation of both (22) and (34) must involve averaging over 
the fluctuating cut- off r\. However, we do not expect this procedure to introduce substantial 
modifications of the obtained results. 

The theory presented here does not involve any multi-fractal assumptions. Still, the quanti- 
tative (numerical) agreements between the two approaches hints on the possibility of some 
qualitative connection. The essentially dynamic theory developed here couples the velocity 
fluctuations at the largest and smallest scales. One may speculate that if a typical structure 
is basically a strongly convoluted sheet with two O(L) linear dimensions and 0(rf) the third 
one, then these structures can loosely be identified with the multi-ifractal sets of the MF 
theory. 
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(a) 




2 4 6 8 

(b) 

Figure 2: a. Computed PDF P(^,r) vs 5 r u/r a (eq (5)) and Log-normal PDF P L {^r,r) 
(eq. (6) ) for the same values of parameters, b. (P — Pl)/Pl vs 5 r u/r a . 
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Figure 3: From Gotoh and Nakano (Ref.[16]). The top curves show the right ("rhs") and 
left ("Ihs") sides of equation (8) for n = 4; 6; 8. When the pressure contribution was added 
to the "rhs", very accurate equality rhs = Ihs has been reached (dotted line on top of the 
"Ihs" curve). The bottom curves are R 2n - 
Fig. 3c is scaled by factor 10. 
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Figure 4: Probability density Q(r)/r)K, Re) vs rj/r/K where t\k is the Kolmogorov scale, a. 
The shape of the PDF Q(rj, Re) for Re w 0.3. b. for Re = 0.1; 0.3. 
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